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We provide a nonparametric method for the computation of in- 
stantaneous multivariate volatility for continuous semi-martingales, 
which is based on Fourier analysis. The co- volatility is reconstructed 
as a stochastic function of time by establishing a connection between 
the Fourier transform of the prices process and the Fourier transform 
of the co-volatility process. A nonparametric estimator is derived 
given a discrete unevenly spaced and asynchronously sampled ob- 
servations of the asset price processes. The asymptotic properties of 
the random estimator are studied: namely, consistency in probability 
uniformly in time and convergence in law to a mixture of Gaussian 
distributions. 

1. Introduction. The volatility is a key parameter in financial economics- 
mathematics and the recent financial econometrics literature has devoted 
much attention to its computation. As a matter of fact, volatility and co- 
volatility computation is crucial in hedging strategies (in the classical Black- 
Scholes environment or in more sophisticated stochastic volatility models), 
in risk management (e.g., VaR methodologies), in forecasting (e.g., ARCH- 
GARCH models) and for optimal portfolio selection. 

Volatility can be computed through parametric or nonparametric meth- 
ods, see, for instance, the review by [4]. In the first case, the expected 
volatility is modeled through a functional form of market or latent variables. 
Nonparametric methods address the computation of the historical volatility 
without assuming a functional form of the volatility. As volatility changes 
over time, its computation through nonparametric methods concentrates on 
a small time window (a day, a week) and high frequency data are employed. 
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In this paper we propose a nonparametric estimation methodology based 
on Fourier analysis which applies to continuous semi-martingales and is 
mainly designed for measuring instantaneous multivariate volatility, exploit- 
ing high frequency observations. 

Suppose the time evolution of the asset price is a Brownian semi-martingale 
of the form 

d P (t) = a(t, W) dW{t) + b(t, W) dt, 

where W is a Brownian motion. Using ltd calculus, the instantaneous volatil- 
ity of the process p is obtained 

(!) Vol W(f ) := a\t) = Jlj. & (M+flzKffi!) . 

Nevertheless it is difficult to use this formula for computing numerically the 
instantaneous volatility. As a matter of fact volatility measurement takes 
place over discrete time intervals and only a single path of the market price is 
effectively observed; therefore, the conditional expectation , with respect 
to the past information jF t , cannot be computed from the observations. As 
a consequence, volatility is mainly computed over discrete time intervals 
(integrated volatility), relying upon the quadratic variation formula. In fact 
a classical result, essentially due to Wiener, states the following formula 
which holds almost surely 

QV|£(p)= f 1 o- 2 {s)ds, 

J to 

where 

QV|£(p):= lim (p(.t + (k + l)2- n )-p(t + k2~ n )) 2 . 

0<fc<(ti-*o)2 n 

In the recent financial-econometrics literature many results exploit the 
quadratic variation formula obtaining nonparametric estimators of inte- 
grated volatility and co- volatilities. The realized volatility — quadratic varia- 
tion estimators have been intensively studied and used for financial-econometrics 
purposes in a series of papers, among them [5, 8, 25, 43]. In order to handle 
irregularly, even randomly spaced data, the quadratic variation theorem was 
extended [2, 9, 28, 37]. 

Nevertheless when instantaneous volatility and co-volatilities are com- 
puted by exploiting high frequency data, according to quadratic variation 
and co-variation formulae, three bottlenecks appear. 

The first one is specific to the multivariate setting: returns are recorded at 
the highest available observation frequency; therefore, they are asynchronous 
across different assets. The realized covariance type estimators require the 
choice of a "synchronization" method, thus these estimators suffer from 
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a downward bias, when the sampling interval is reduced (known as Epps 
effect, by [16]). Reference [42] provides an analytical study of the realized 
co-volatility estimator in a general framework which includes asynchronous 
trading; recently [20] have proposed an estimator of integrated co-volatility 
which is consistent under asynchronous observations. 

Second, instantaneous volatility computation (1) involves a sort of nu- 
merical derivative. Therefore, the authors, who address the problem of es- 
timating the instantaneous volatility, use a double asymptotic in order to 
perform both the numerical derivative involved in formula (1) and the ap- 
proximation procedure (see [14, 18, 37]). However, the empirical derivative 
required to get the instantaneous volatility gives rise to strong numerical 
instabilities; moreover, none of these authors address the problem of asyn- 
chronous observations in the estimation of instantaneous co-volatility. A 
different approach through wavelets' series is developed in [19] for deter- 
ministic volatility function. Finally, functional methods, which are local in 
space, have been developed for estimating the volatility as function of the 
underlying state variable level; see [1, 7, 17, 26]. 

Third, the limiting argument n — > oo relying on high frequency data can- 
not be effectively accurate due to market microstructure effects (such as 
discreteness of prices, bid/ask bounce, etc.), which cause the discrepancy be- 
tween asset pricing theory based on semi-martingales and the data at very 
fine intervals. The impact of market microstructure noise on the realized 
volatility estimator has been studied, among others, in [3, 15, 41, 44], and 
modifications of realized variance estimator have been proposed to correct 
the bias due to microstructure noise, for example, the two-scales estimator 
by [43] and the realized kernels by [10]. 

The estimation procedure proposed in this paper is immune from the 
first and second bottlenecks highlighted above, in virtue of its own defi- 
nition. Fourier estimator uses all the available observations and avoids any 
"synchronization" of the original data, because it is based on the integration 
of the time series of returns rather than on its differentiation. As well as re- 
garding the third bottleneck, which concerns the robustness of the estimation 
procedure under microstructure noise, the Fourier estimator is competitive, 
in terms of mean squared error, with the methods specifically designed to 
handle market microstructure noise with only a slightly higher bias; see [36]. 

The first step in the Fourier analysis methodology for the computation 
of co- volatility is [32] , where the instantaneous multivariate volatility func- 
tion is expanded into trigonometric polynomials, whose coefficients depend 
on the log-return processes. The efficiency of the Fourier series procedure 
for computing and forecasting the integrated volatility has been analyzed in 
[12, 13, 21, 29, 38]. Moreover, an application of Fourier method to the com- 
putation of correlation dynamics of futures markets is performed in [24]. 
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In this paper we extend the results of [32]: we prove a general identity 
which relates the Fourier transform of the co-volatility function with the 
Fourier transform of the log-returns under the hypothesis that the volatility 
process is square integrable. From this identity we derive an estimator of 
instantaneous co-volatility based on a discrete, unevenly spaced and asyn- 
chronously sampled asset prices. Finally, we prove the consistency in prob- 
ability uniformly in time of the proposed estimator and the convergence 
in law to a mixture of Gaussian distribution. As a byproduct, the inte- 
grated volatility and co-volatilities estimators are obtained by considering 
the k = coefficient of the Fourier expansion. Our theory implies that the 
Fourier estimator of integrated co- volatility is consistent under asynchronous 
observations in the absence of microstructure noise. 

It is worth noting that the Fourier estimator has two adjustment param- 
eters: the number of coefficients computed for the price and the number of 
Fourier coefficients used for the reconstruction of the volatility. Balancing 
the choice of these two parameters according specific situations, the estima- 
tor can act as an optical microscope with different magnifications: in the 
case where microstructure is supposed to be produced by irrelevant noise, 
low magnification will be used. In both cases, the whole set of observed data 
are employed. In summary, the Fourier estimator is designed specifically for 
high frequency data: by cutting the highest frequencies, it uses as much as 
possible of the sample path without being more sensitive to market frictions; 
see also [36, 38]. 

Finally we stress the following point: the potentiality of Fourier estimation 
method relies on the fact that it allows to reconstruct the volatility as a 
stochastic function of time in the univariate and multivariate case. This 
property makes possible to iterate the volatility functor and to compute the 
volatility of the volatility; see [33] for an early recognition on this point. In 
other words, we can handle the volatility function as an observable variable. 
For instance, this feature is essential when a stochastic derivation of volatility 
along the time evolution is performed as in contingent claim pricing-hedging 
(see [11, 35]), or when we study the geometry of the Heath-J arrow-Morton 
interest rates dynamics given the observation of a single market trajectory 
[34]. 

The paper is organized as follows. Section 2 contains the main result: we 
derive the Fourier transform of the volatility function as the Bohr convolu- 
tion product of the Fourier transform of log-returns. In Section 3, given a 
discrete, unevenly spaced and asynchronously sampled observations of the 
price processes, we prove the consistency uniformly in time of the random 
estimator of the co- volatility function by considering a fixed time span and 
increasing the frequency of the observations. In Section 4, the weak conver- 
gence in distribution of the co-volatility process to a mixture of Gaussian 
laws is proved. In Section 5 we analyze the implications of our results for 
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the Fourier estimator of integrated volatility and co-volatility, even in the 
presence of microstructure noise effects. Section 6 concludes. Proofs are in 
Section 7. 

2. Instantaneous co- volatility computation. Let p(t) = ■ ■ ■ ,p n (t)) 

be a continuous semi-martingale satisfying the following Ito stochastic dif- 
ferential equations 

d 

(2) d I j(t) = Y,<4(t)dW i + V(t)dt, j = l,...,n, 

i=i 

where W = (W 1 , . . . , W ) are independent Brownian motions on a filtered 
probability space satisfying the usual conditions, and cr* and b* are adapted 
random processes satisfying 

\2 



(H) 



E 



E 



(b l (t)Ydt 



{oi{t)fdt 



< oo, 



< oo, 



1, 



,d, j = l, 



,n. 



In financial applications the process p{t) represents the log-price of assets 
and from the representation (2) we define the volatility matrix, which in our 
hypothesis depends upon time: 

(3) ^w=E°!(^w- 

8=1 

For notational simplicity we will refer to the case of two assets whose 
prices are (p 1 (t) , p 2 (t)) . We compute the so called instantaneous volatilities 
E"(t) for j = l,2 and co-volatilities S 12 (t) = S 21 (t). 

By change of the origin of time and rescaling the unit of time we can 
always reduce ourselves to the case where the time window [0, T] becomes 
[0,27r]. We recall some definitions from harmonic analysis theory (see, e.g., 
[30]): given a function (j) on the circle S , we consider its Fourier transform, 
which is defined on the group of integers Z by the formula 

F(<j>)(k) := — <t>($) exp(-ifa?) d$ for k € Z. 



We define 



F(d(j))(k) := — ( exp(-i/h?)#(tf), 

2-7T 7l0,27r[ 



and, using integration by parts, we remark that 



1 

2^ 



(0(27r)-0(O))-JP(#)(fc) 



G 
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Given two functions <!>, on the integers Z, we say that the Bohr convolution 
product exists if the following limit exists for all integers k 



1 N 

(4) (** B *)(*):= lim — — HsMk-s). 

TV— >oo ZiV + 1 * — ' 

s=—N 

The subindex B is the initial of Bohr, who developed similar ideas in the 
context of the theory of almost periodic functions. 

The main result of this section is formula (5) on which the construction 
of the volatility estimator relies. Theorem 2.1 contains the identity relating 
the Fourier transform of the price process pit) to the Fourier transform of 
the volatility matrix 

Theorem 2.1. Consider a process p satisfying the assumption (H). 
Then we have for i,j = 1, 2: 

(5) ±-F(XV)=F(d P i )* B T(dpP). 

The convergence of the convolution product (5) is attained in probability. 

The proof is based on the ltd energy identity for stochastic integrals. A 
preliminary step proves that the drift component b* (t) of the semi-martingale 
p(t) gives no contribution to the formula (5). The latter derives from the 
following result. 

Lemma 2.2. Consider a function u defined on Z and define the Bohr 
L? -seminorm 



1 N 

u\\ BL 2 :=limsup ( u ( k )) 2 ■ 



Then 

(6) \\u *b v\\l°o < \\u\\ BL 2\\v\\ BL 2. 
In particular, it follows that 

(7) || u IIsl 2= =^ (u*bv) = for all u. 



By Theorem 2.1 we gather all the Fourier coefficients of the volatility ma- 
trix by means of the Fourier transform of the log-returns. Then the recon- 
struction of the cross- volatility functions ^(t) from its Fourier coefficients, 
can be obtained as follows: define for i,j = 1,2 

$jv(fc) : = ^{dp l ){k) for|/c|<2iV and otherwise, 
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and, for any \k\ < N, 

==«iirn:E*wW*ir(*--)- 

sGZ 

If E lJ (t) is continuous, then the Fourier-Fejer summation gives almost ev- 
erywhere 

(8) S«(t)= Urn £ (l-^l)^)exp(i^) for all t e (0, 2vr). 

Remark 2.3. We note that in the case of the volatility functions T, ll (t), 
the approximating trigonometric functions appearing in the sequence (8) 
are positive. In fact, by the positivity of the Fourier transform of the Fejer 
kernel it is enough to show the same result for the function 

QNit) :=J2¥ N (k)exTp(ikt). 
kez 

Then we have 

(2N + l)Q N (t) = J2 exp(iH) £ ^ l N (sW N (k - s) 
fcez sez 

and this is equal to 

2 



]T eMKs-s')tW N (s)& N (s') 



^exp(ist)$V(^ 
sez 



using the fact that s) = & l N (s), where denotes the conjugate func- 

tion of <I>\ 

In the case of a single asset, the computation of the volatility function 
follows by Theorem 2.1. For completeness we state the result. 

Corollary 2.4. Suppose that the log-price process p(t) follows the ltd 
process 

dp{t) = b{t)dt + a{t)dW{t), 

where W is a standard Brownian motion and a, b are adapted random pro- 
cesses satisfying assumption (H). Then the Fourier coefficients of the volatil- 
ity function o- 2 (t) are obtained, for any integer k, by 

±f(a 2 )(k)=lim— 3— ]T F{dp){s)F{dp){k-s), 

\s\<N 

where the limit is attained in probability. 

The instantaneous volatility function o~ 2 (t) can be obtained analogously 
with (8), using the Fourier-Fejer summation. 
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3. Consistency of the Fourier estimator of instantaneous co- volatility. In 

this section we prove that the estimation methodology proposed does not 
induce any synchronization bias and the Fourier estimator of instantaneous 
co- volatility is consistent. As a particular case the consistency of the inte- 
grated co-volatility estimator follows. 

Let S ni ■= {t} m ,i = 1, . • • ,ni} and T n2 := {t 2 n2 ,j = 1, ... ,712} be the trades 
for the asset 1 and 2, respectively. For the ease of notation we often sup- 
press the second index ni,n,2- For simplicity suppose that both assets trade 
at to = an d i ni = t„ 2 = 271". It is not restrictive to suppose that ri\ = ri2 =: 
n. Denote for k = 1,2, p k (n) := max </i< n _i \t\, 1 — t^\ and suppose that 
p(n) := /^(n) V p 2 (n) — > as n — > 00. 

Consider the following interpolation formula 

n-1 



pKt) :=£^(4)% iU (*)- 



i=0 

Analogously for Pn(t). We will indicate // := [tj,t} +1 [ and J| := [iy,t? + i[ and 
the returns by S I i(p 1 ) :=p 1 (^ +1 ) -p 1 ^}) and Sj2(p 2 ) :=p 2 (t 2 +l ) -p 2 (t 2 ). 
For any integer k, \k\< 2N, let 



(9) 



^ n— 1 

(dpi) ■= — EexpH^ 1 )^ 1 : 



2 n— 1 

Ck{dp 2 n ) -=^Y1 ex P(- ikt2 j) 6 Jj(P 2 ) 

j=0 



moreover define 



1 



2tt 



(10) c k {T} 2 ):=— e- ikt E 12 (t)dt. 

ZTT Jo 

Let for any \k\ < N 

2rr 

(11) a k (N,pl,pl) := — — - 2 Csidp^cksidpl). 

\s\<N 

Finally, define 

(12) S^(t) := E f 1 " §W^Pn,Pn)e ifci . 

|fc|<AT V 7 

The random function T,^ N (t) will be called the Fourier estimator of the 
instantaneous co- volatility £ 12 (i). In the following theorems we prove the 
convergence in probability of the random function (12) to the co-volatility 
function S 12 (t) uniformly in t, as N and n go to infinity. 
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Before stating the result, some comments are needed. We note that the 
effective reconstruction of X 12 (i) is realized as superposition of three limits. 
The first limiting procedure consists in the convergence as ./V goes to infinity 
of the Fejer sum 

E (l-§)exp(i^)c fe (S 12 ). 

\k\<N K J 

To this aim we remark that, given a continuous function 
its modulus of continuity defined as 

uv(A):= sup 

|0-0'|<A 



denote by u^(A) 



and let Ck(4>) '■= J~(cft)(k); then it holds 

|fc| 



(13) 



sup 

t 



E 

\k\<N 



1 



N 



ex.p(ikt)ck((j)) 



< 



^0 



(see [40] for the converse of this inequality). Therefore, if we assume that 
some a priori information on the modulus of continuity of S(t) is available, 
from (13) we see that the determinacy of £ 12 (t) in the L°°-norm can be 
obtained from the knowledge of a finite number of its Fourier coefficients. 
In the sequel suppose that S*- 7 (t) satisfies 

(Hi) esssup US 1 - 7 ||ioo < oo, where US 1 - 7 ||^oo := sup |£ lJ (t)|. 

t 

The determination of Cfc(E 12 ) depends upon the superposition of two lim- 
iting procedures, letting n go to infinity. First, we approximate c^dp 1 ) by 
c k{dPn)- The second limit procedure appears in the definition of the Bohr 
convolution. Next, two lemmas provide the required inequalities to obtain 
the consistency result. 

Lemma 3.1. For any integers k,s, with \k\ < N and \s\ < N, and any 
n, the following inequality holds: 

E [\cs(dpl)c k _ s (dpl) - Csidp^Cksidp 2 )] 2 ] 

(14) 

< CiV 2 p(n) 2 esssup HS 11 )^^ esssup ||S 22 ||x,oo, 
where C is a suitable constant. 



LEMMA 3.2. For any k, \k\ < 2N , the following inequality holds: 
E 



(15) 



2iV + 



< 



- E Csidp^Cksidp 2 ) -c fc (S 12 ) 
|s|<JV 



ll^llll ll^22|i 

esssup ||2j 1 1 l°° esssup || h \\lo 



2N + 1 
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Theorem 3.3. Given Cfc(<ip„) and Ck{dp^) defined in (9) and c k (T, 12 ) 
defined in (10), let for any \k\ < N 

(16) a k (N,p n ,pl) := —— ^ c s (dp n )c k _ s (dp n ). 

\s\<N 

Suppose that T,(t) is continuous and Np(n) — > as N,n—> oo. Then, for any 
k, the following convergence in probability holds: 

(17) lim a k (N,p 1 n ,p 2 n ) = c k (T, 12 ). 

n,N— >oo 

Finally, we conclude by [40]: 

Theorem 3.4. Let T^ N {t) defined in (12). If E(f) is continuous and 
N p(n) — > 0, the following convergence holds in probability 

lim sup |S^(t)-S 12 (t)|=0. 

n,A r ->oo <t<27r 

In the univariate model specified in Corollary 2.4 the Fourier estimator 
of instantaneous volatility can be obtained as follows. 

Consider a discrete unevenly spaced sampling of the process p. Fix a 
sequence S n of finite subsets of [0,2-7r], let S n := {0 = to,n < h,n < •• < 
*fc n ,n = 27r} for any n > 1 such that /?(n) := maxcK/^/^-i - th,n\ -> 

asn^oo. Denote, for any j, <5j(p) :=p(t J+ i jn ) -p(t jtn ). 

Define, for < JV, 

27T 

(18) a k (N,p n ) '■= 2N + 1 Yl c s(dPn)ck-s(dp n ), 

\s\<N 

where for any integer k, \k\ < 2N , c k (dp n ) are the Fourier coefficients of the 
log-return process, that is, 

Ck{dpn) = — Y] exp(-iktj >n )6j{p). 

Corollary 3.5. Consider for any n,N the random function of time 

(19) Zl, N (t)--= E U-^)MN,Pn)eM^t), 

\k\<N K J 

where a k (N,p n ) is defined in (18). Suppose cr 2 (t) is continuous. 

(i) For any integer k the following convergence in probability holds 
lim a k (N,p n ) = c k (o 2 ). 

n,N—*oo 
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(ii) It holds in probability 

lim sup \B 2 N {t) - a 2 {t)\ = 0. 

n ' N ^°°t€lO,2n} 

Remark 3.6. In the construction of the Fourier estimator we can con- 
sider observations {i^} which are random but independent of the observed 
processes p k , k = 1,2. In fact by splitting the probability space, this kind of 
sampling reduces to a deterministic sampling. Therefore all our study has 
been made in this context. 

4. Limiting distribution. In this section we prove the weak convergence 
of the instantaneous co-volatility estimator defined in (12) to the process 
£ 12 (i). As a corollary the asymptotic normality for the error of the integrated 
co- volatility estimator follows. The asymptotic results are obtained as the 
mesh of the partition p(n) — > and the interval [0, 2ir\ remains fixed. Let 
N be the number of Fourier coefficients considered in the reconstruction of 
and suppose that N — > oo. 

We first consider the univariate model specified in Corollary 2.4 and the 
estimator of the volatility function N (t) given by (19). As we suppose to 
have an irregularly spaced grid, some conditions are needed. The following 
assumptions from [37] are sufficient: 

(A) (i) p(n) -> and k n p(n) = 0(1) , 

(") H n (t) := ^> +1 ' n % /kn >H(t) as oo, 

(iii) H{t) is continuously differentiable. 

In the sequel we will denote by tp the conjugate of the complex function 

Theorem 4.1. Let o~ n ]\[{t) defined in (19). Assume that as n — > oo then 
p(n) — ► 0, — ► oo, and assumption (A) is satisfied. Then, for any function 
h G Lip(a), a > \, with compact support in (0, 2ir), 

(p(n)r 1/2 f 2 \{t)(^{t)-a 2 (t))dt 
Jo 

converges in law to a mixture of Gaussian distribution with random variance 
2f 27T H'(t)h 2 (t)a A {t)dt, provided that p(n)N 2a ^oo. 

The result in Theorem 4.1 derives from a central limit result of [27], gen- 
eralized for unevenly spaced data in [37] (Proposition 1), after the following 
representation theorem holds. 
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Theorem 4.2. Under the hypothesis of Theorem 4-1 the following result 
holds 

/■27T 

(pin))- 1 / 2 h(t)(al N (t)-a 2 (t))dt 
Jo 

kn — 1 rtj+1 n f ft \ 

(20) = (p(n))" 1/2 2 J2 M*i,») / / a( S )dW(s) W(t)dW(t) 

+ o p (l), 

where the function h^ is defined as hjsr(t) := J2\k\<N(^ ~ ^N) c k(h)e lkt and 
o p (l) is a term which goes to zero in probability. 

Remark 4.3. The above study of the asymptotic properties of the esti- 
mator also provides a way to obtain the order of magnitude of the number 
of Fourier coefficients, N, to be included in the estimation given the number 
of observations, n. The asymptotic properties require that p(n)N — > and 
p(n)N 2a — > oo, if p(n) is the mesh of the partition. Similar conditions appear 
in order to derive the limiting properties of Bartlett type kernel estimator 
(see [43]) for the relative growth of the number of auto-covariances with 
respect to the number of data. 

The asymptotic result in the bivariate setting needs to take into account 
for the asynchronicity of the returns. Then assumption (A)(ii) must be mod- 
ified according to 

(A) (if) H n (t) := £ Ztj +1 At* +1 <Sl +l A t 2 +l - tj V ^) 2 / {t i vt|<t i +iAt?+i} - 
H(t) as n — > oo. 

Theorem 4.4. Let T, 12 N (t) defined in (12). Suppose p(n) := p l (n) V 
p 2 (n) — > 0, N — > oo, and assumption (A) (i), (ii'), (hi) hold, then, for any 
function h G Lip(a), a > |, with compact support in (0,27r), 

(21) ( P (n))- 1 ' 2 h(t)(^ 2 N {t) - S 12 (t)) dt 

Jo 

converges in law to a mixture of Gaussian distribution with variance Jq W H'(t) x 
/i 2 (t)(S n (t)S 22 (t) + (S 12 (t)) 2 ) dt, provided that p(n)N 2a -> oo ; if a > §, and 
p(n)N 4 / 3 -f 0. 

In order to apply central limit result of [27] (see also Section 4 in [9]), the 
following representation theorem is required. 



A FOURIER TRANSFORM METHOD FOR NONPARAMETRIC ESTIMATION 13 

Theorem 4.5. Under the hypothesis of Theorem 4-4> the following result 
holds 



o p (l) is a term which goes to zero in probability. 

5. Fourier estimator of integrated volatilities in the presence of microstruc- 
ture noise. The financial econometric's literature mainly focus on the inte- 
grated volatility and co- volatility, that is, Jq S 1 - 7 (t) dt, for i,j = 1,2, where 
[0, T] is a fixed time horizon, for example, a day. In fact the computation of 
spot volatility by using quadratic variation formula is quite instable, because 
it involves a numerical derivative. 

In the context of Fourier estimation methodology the integrated volatility 
and co- volatility are computed by considering the k = coefficient, respec- 
tively, in the formulae (11) and (18). The theory presented in this paper 
is stronger because we recover pathwise the entire volatility curve from the 
observation of an asset price trajectory. Nonetheless, due to the importance 
of the integrated (co-) volatility for financial econometrics purposes, in this 
section we explicitly state the results concerning integrated volatilities; in 
particular, the properties of Fourier estimator in comparison with different 
integrated volatility estimators proposed in the literature can be analyzed 
in this context. 

Consider first the univariate case. Suppose that the price process p(t) is 
a continuous semi-martingale satisfying the hypothesis stated in Corollary 
2.4. Observe that by definition 



Jo 

therefore the integrated volatility over a given period (always reducible to 
[0,27r]) is computed by Corollary 2.4 
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More precisely, using the same notation as in Section 3, from (22) the 
Fourier estimator of integrated volatility over [0, 2tt] is defined as follows 

(23) ^' Ar:= 2^+T c s (dp n )c^ s (dp n ). 

\s\<N 

The consistency in probability uniformly in time of the Fourier estimator 
of the integrated volatility follows immediately by Corollary 3.5. The effi- 
ciency of the Fourier procedure to obtain estimator of integrated volatility 
is studied in [13, 21]. We remark that the Fourier estimator (23) can be 
expressed as 

(24) a 2 n>N = ££ D N (tj - t' J )5,(p)5 f (p), 

j=ij'=i 

where -Djv(i) denotes the rescaled Dirichlet kernel defined by 

(25) DN (t) : =— l_ v «- = 1 ^m±iim. 

V ' W 2iV + l,f-'„ 2N + 1 sin (t 2) 

\s\<N v ' ' 

From (24) different features of the Fourier estimation method are high- 
lighted. First, the Fourier estimator uses all available data by integration, 
thus incorporating not only the squared log-returns but also the products of 
disjoint log-returns along the time window. Second, the convolution product 
weights the auto-covariances at any given frequency. In fact it is easily seen 
from (24) that 

o* tN = RV+Y, D N (tj - gSj^Sfip), 

where RV denotes the realized volatility estimator, namely X)j=i(^?'(p)) 2 ■ 
The contribution of various order auto-covariances has early been considered 
by [44] and recently used to correct the bias of the realized variance type esti- 
mators in the presence of microstructure noise, in particular the subsampled 
estimator by [43] and the realized (subsampled) kernels by [10]. Different 
weighting functions are also considered in [6, 29] who express many volatil- 
ity estimators proposed in the literature in an unified way as a quadratic 
form. In the analysis of [23, 29] a peculiar feature of Fourier estimator is not 
efficiently addressed: the Dirichlet kernel appearing in the Fourier volatil- 
ity estimator depends on the number of frequencies N, besides the number 
of auto-covariances and the time lag between two observations. The im- 
portance of the choice of the frequencies has been empirically noticed by 
[38]; the authors make an empirical comparison, through Monte Carlo sim- 
ulations and using high frequency market data, between some estimators 
of integrated volatility. Their analysis shows that the Fourier estimator is 
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superior to the realized volatility and the wavelet estimators and, even com- 
pared to some bias correcting methods for microstructure noise, namely the 
kernel-based estimator by [44] and the related unbiased estimator proposed 
by [22], the Fourier estimator provides smaller root mean squared error while 
having only slightly higher bias. Recently [36] prove that, while the MSE 
of the realized volatility estimator under microstructure noise diverges as 
the number n of observations increases, the MSE of the Fourier estimator is 
substantially unaffected by the presence of microstructure noise by choosing 
in an appropriate way the number of Fourier coefficients to be included in 
the estimation, as indicated explicitly by the MSE computation. The op- 
timal MSE-based Fourier estimator turns out to be very efficient, even in 
comparison with methods specifically designed to handle market microstruc- 
ture contaminations; more specifically, the Fourier estimator is competitive 
in terms of MSE for high sampling frequencies up to 30 seconds. We re- 
mark that in the presence of microstructure noise it is useful to consider the 
following more stable version of Fourier estimator 

(2ir) 2 N / Is! N 
(26) ^iV-Trrr E (l-^)cs(dpn)c-s(d Pn ). 

s=-N ^ ' 

As for the multivariate case, the integrated co-volatility is recovered as a 
byproduct of formula (5) through the following identity 

2vr^(S 12 )(0)= f 2T J: 12 (t)dt. 
Jo 

Therefore the Fourier estimator of integrated co-volatility is defined as 

1 n— In— 1 

* + \s\<Ni=0j=0 

The consistency of the Fourier estimator of integrated co-volatility with 
asynchronous observations for different assets follows from Theorem 3.3. 
Moreover by (14) and (15) we get an estimate of the mean squared error 
(MSE) of the Fourier integrated co- volatility estimator 

MSE 9l2 <(c 1 - + C 2 N 2 p(n) 2 ) ess sup ||£ n ||£oo ess sup ||£ 22 ||£oc, 

n,iV \ N J 

C\ and C2 being constants. By choosing N such that = N 2 p(n) 2 (up to 
a constant) we obtain that N = p(ra) _2//3 . Therefore the mean squared error 
of Fourier estimator of co-volatility is of order 0(p(n) 2 / 3 ). 

All the realized covariance type estimators rely on "synchronization" 
methods; therefore, they suffer from a downward bias when the sampling 
interval is reduced (known as Epps effect by [16]). Recently [20] have pro- 
posed an estimator of integrated co-volatility which is consistent under 
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asynchronous observations. This estimator is defined by adding the cross- 
products of all fully and partially overlapping returns. [39] prove that the 
estimator of [20] is biased in the presence of microstructure noise and its vari- 
ance diverges as the number of observations increases. We note that Fourier 
estimator of integrated co-volatility defined in (27) has the advantage that 
from one side it is consistent in the absence of microstructure contamina- 
tions, but with the optimal choice of iV (as explained in the univariate case), 
it is more robust in the presence of microstructure noise effects. 

In summary, the Fourier estimator exploits more information in the data 
without being affected by a severe bias coming from microstructure contam- 
ination, simply by choosing the lower frequencies. 

6. Conclusions. In this paper we have proposed a nonparametric estima- 
tor of the instantaneous volatility and co-volatility function in the context 
of continuous semi-martingale models, given the knowledge of the Fourier 
transform of the log-returns. The method applies both to univariate and mul- 
tivariate setting and allows to consider unevenly spaced and asynchronous 
observations. Moreover we have studied the statistical properties of the 
Fourier estimator by providing consistency and asymptotic normality. We 
emphasize two main features of the proposed estimation methodology: the 
first one is that the definition of Fourier estimator is designed specifically to 
handle the computation of co- volatilities without any manipulation of data, 
due to the fact that the Fourier method requires an integration of all data 
at disposal for any asset; the second one is that, with the availability of high 
frequency data, we are able to reconstruct the co-volatility as a stochastic 
function in an effective way, which permits to handle the volatility function 
as an observable variable in the econometric and financial applications. 

7. Proofs. In the sequel we will denote by ip the conjugate of the complex 
function ifj. 

Pro of of Lemma 2.2. For a € Z denote u a (k) := u(k + a) and u(k) := 
u(—k). Observe that 



(u *b v)(a) = (u\v 
where {u\v)bl 2 denotes the scalar product 



a 



) 



BL 2 i 



1 



N 



{u\v) BL 2 := lim 



iV-»oo 2N + 1 



k=-N 



By Cauchy-Schwarz inequality it holds: 



1 N ( 1 N 

y u(k)v(k) < — - — y ( u (k)) 




1/2 
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This implies that for all a 

(u* B v)(a) < \\u\\ BL 2\\v" 
Because ||#~ a || S £ 2 = |H|bl 2 ) the result is proved. □ 



- — ail 

I RE- 



PROOF OF Theorem 2.1. Without loss of generality we suppose that 
b* = 0. In fact, consider the process p m = {Pm,Pm) given by dp> m {t) = 
Ya=i °\ (*) dW l (*) : Pm has the same volatility function than the process p. 
We have that for % = 1, 2: 

Hdp i ){k)=Hdp i m ){k)+HV){k). 

For notational simplicity we omit now the upper index j = 1,2. With the 
notation $ m (fc) := J : (dp m )(k), <&b{k) := T(b)(k), we now prove that 

(* m + $ 6 ) *B ($m + *&) = $m *B $m, 

because the terms $6 $>{,, $6 $ m and <E> m *_b give zero contribution. 
First, we consider &b*B We have 

1 " 

($&* B $ b )(fc)= lim 51 *6(a)*6(fc-«)- 

s=— TV 

It follows by the Lemma 2.2 that: 

(28) ||$b*B *6||l« < ||*6||Bjy>l|*6||Br2- 
Moreover by Plancherel theorem it holds: 

+°° /•27T 1 

(29) £ |^(6)(fc)| 2 = y -|^)| 2 dt. 

fc=— oo 

Therefore by (29) it follows: 

N 



(30) ||$ 6 ||| L2 = Jim J] |^(6)(A:)| 2 = 0. 



The fact that the terms &b*B &m an d $ m *B &b gi ye zer ° contribution 
follows by the Lemma 2.2 and (30). 

Under the hypothesis that b* = 0, then p is a martingale. Introduce the 
complex martingales for any integer k and i = 1,2 



r*(*) :=7^ fzM-iks)dp\s) 
2ir Jo 



Then by definition r^.(27r) = J r (dp l )(k). Compute by Ito formula the stochas- 
tic differential 

d(riTi)(t) = Z ij (t) exp(-i(jfc + r)t) dt + T{(t) dTi(t) + T*(t) <ir J fc (t), 



18 

we obtain 

(31) 

where 



o 

Fix an integer N > 1 and define for any integer q with \q\ < N, 

N 
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r fc (27r)P r (2vr) = ^?(P 3 )(k + r) + BP(k, r), 



W(k, r) := I ^ Ti(t) dTi(t) + K(t) dTUt). 



(32) 



1 



< j ( N ) = WTl £ n +s (2^)ri s (2vr). 



s=-N 



By definition (4), it results that, for any integer q 

lim 1 ^{N) = (T(dp i )* B T(dp>))(q). 



Moreover by (31), we have: 



(33) 
where 



1 



N 



q "' / 2vr 
r 2w N 



N> 



SvTT /„ E rj-i(*X(*) + rf 

T JU l=—N 



(t)dr*_,(t). 



We prove that i?^y converges to in probability as iV — ► oo. By symmetry 
we are reduced to study 

i N r 2w . r t 2 



Introducing the rescaled Dirichlet kernel 



1 



A' 



D N(t) := 2N+1 ex P( ist ) 



s=-JV 



1 sin(iV + l/2)t 
2N + 1 sin(i/2) ' 



we have that 



(34) A 
Define now 



iV 



47T 2 Jo 



2 IT 



dp>{t 2 



exp(iqti)D N (ti - t 2 ) dp l (t\ 



a(t 2 ) := f 2 coB(qti)D N (t! - t 2 ) V 4(h) dW k (h 
Jo f~t 



k=l 
2 



/3(t 2 ) := f * sm(qh)D N (h - t 2 ) V 4(h) dW k (h 
Jo fc=i 
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By ltd energy identity for real stochastic integrals we have (e.g., [31]) 



167T 4 ^[|^| 2 ] = V^[/ 27r ( 
k=l U ° 

By Cauchy-Schwarz inequality we have 
(167r 4 E[[^| 2 ]) 2 

2 r r 2TT 



a 2 (t 2 ) + (5 2 (t 2 ))(a{(t 2 )fdt 2 



<4J2E 



k=l 



K(t2)) 4 dt 2 



x < E 



£ n ^ 2 D N (h -t 2 )co S (qt 1 )Y J 4(ti)dW k (t 1 ) ) j dt 2 



+ E 



( f 2 D N (h - t 2 ) sm(qti) V 4(h) dW k (h 
Jo V° k=i 



dU 



Because (cr , k (t)) 4: dt] < oo by hypothesis, we evaluate the other term 

using Burkholder-Gundy inequality 

4 

E 



jf 7r ^ 2 Djv(ii-t2)cos(^i)^4(ti)^ fe (ti)J dt 2 



<AE 



2tt r t 2 







[ 2 D%(t 1 -t 2 )Y l (<(h)) A dt 1 dt 2 ; 



then making a change of variable t\ = u, t 2 —t\= v, we get 

- r2n 2 1 /-27T 

= 8E / 5>£(«)) 4 du 
Uo k=i 

We remark that |Z)jv(^)| < 1> therefore 



o 



D%(y)dv. 



2tt 



2jt 



D%(v)dv< / D 2 N (v)dv 



2tt 



2iV + l 

The last equality is a consequence of Plancherel equality. Therefore letting 
N — > oo from (33) we get the convergence in probability. □ 

PROOF of Lemma 3.1. For any i = 1,2 it holds 

1 /" 27r 



c k (dp l n ) - c k {dp % 



2vr Jo 



P(t)d P l (t), 



where 



n-l 



0(t) :=^e ifet ?(l- e ifc (*-*!))I^^ +i[ (t) 
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Now by Burkholder-Gundy inequality we have 
(35) EWckidpi) - c k {d P l )\ 4 \ ^CfcV^esssupll^Hloo, 

where C is a suitable constant. Therefore combining (35) and the inequality 
E [\ c s{dpi)c k - S (dpl) - c s {dp l )c k ^ s {dp 2 )\ 2 ] 

<2(i?[| Cs (dp 1 )| 4 ] 1/2 ^[|c fc _ s (^ ) _ Cfc _ s(dp 2 )|4] l/2 

+ ^[| Cfc _ s (dp 2 )| 4 ] 1 /2 i i;[| Cs (^)- Cs (V)| 4 ] 1/2 ), 
we complete the proof. □ 



Proof of Lemma 3.2. By Ito formula 
2tt 



E 



2N + 



- Csidp^c^sidp 2 ) - c fc (S 



12- 



\s\<N 



<2\E 



2n 



+ E 



o 

1 f 2n 
2^ 



j^DNih-t^dp'ih 



dp 2 (t 2 ) f 2 e ikt 'D N (t 2 - h) dp\U 
Jo Jo 



2i 



where -Djv(i) is the rescaled Dirichlet kernel. Let 



A N (t)= / e^D N (t-v)dp 2 (v), 
Jo 

then by Ito energy identity and Holder inequality 



2tt 



E 



as j^D 2 N (v)dv = J ^ ri . □ 



dp\t 2 )A N (t 2 ) 



E 



2tt 



A 2 N {t 2 )Y} l (t 2 )dt 2 



2tT ii iiii II o o ii 

< ess sup £ /.ooesssup S 



Proof of Theorem 3.3. Using the inequalities obtained in the Lem- 
mas 3.1 and 3.2 it follows 

<12\ 



E 



2tt 



,7 E c s (dpi)c fc _ s (dp 2 ) -c fe (S 

|s|<JV 
1 



< \Ci- 



■2N + 1 

where Ci and C2 are constants. □ 



+ C 2 N p{n) esssup||S ||x,oo esssup ||£ ||z,°°, 
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Proof of Theorem 3.4. In the previous theorem we have proved for 
any fixed k the convergence in probability of each coefficient a k {N ^^^p 2 ^ 
to the Fourier coefficient Cfc(S 12 ), as n and N go to oo. Therefore the thesis 
follows by applying (13). □ 

Proof of Theorem 4.2. By using the Fourier series decomposition 
we have 

2tt 



M 1 

|fe|<iV v 



N 



c k (h)a k (N,p n ) - E c k (h)c k (a 2 ) - E c k (h)c k (a 2 ), 

\k\<N \k\>N 



where c k {f) denotes the fcth Fourier coefficient of the function /, a k {-) is 
defined in (18). By Cauchy-Schwarz inequality and applying (13) we have 



(36) {pin))" 1 ' 2 £ {l-^) Ck { h n{a 2 )<C{p{n)r l / 2 N~^a 2 \\ L ,, 
\k\>N^ 



where C is a suitable constant. As p(n)N 2a — > oo, then (36) goes to as n, 
N ^oo. 

Now we consider 

\k\ 

)c k (h)aj:(N,p r 



2- E i 

\k\<N V 



N 



h n 1 



^h N {t hn ){5 J {p)) 2 + 2 E (l- l -^)c k (h)A N (k), 

j=0 \k\<N K J 



where 
(37) 

and 

It holds 



h N (t):= E 



ikt 



\k\<N 



A N(k) := E ^'"^ivTl ^ ^V,«-^)^(p) 5 .,(p). 



3<3' 



\h\<N 



E 



E (l- l -§)c k (h)A N (k) 



(38) 



|fc|<JV 



E h N(tj,n)DN(tj',n ~ tj,n)5j(jp)5ji (p) 
3<f 
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2tt 



< C\ ess sup ||<t 2 || 2 



2JV + 1' 



where C\ is a constant. It follows that (p(n)) 1//2 .E[| X)|fc|<jv(l ~~ w) c k(h) 
A]v(k)\ 2 ] goes to in probability as n,N — > oo. Therefore it holds 

(p(n))- 1 /2 2 vr £ (\-M) Cfc(/l) ^ (i V )Pn ) 



(39) 



|fc|<JV 

/CjT, 1 



(p(n))- 1 / 2 £ ^fen)(^(p)) 2 + 0p (l) 
J'=0 



Finally it remains to consider 

h n 1 



3=0 

It is enough to study 



E h N(tj,n)(Sj(p)) 2 - Ck{h)c k {a 2 ). 

\k\<N 



k n 1 
3=0 



2tt 



h N {t)a 2 {t)dt 



by It 6 calculus this is equal to 

fen— 1 / r t 



2 J2 M*j>) ( 1 3+1 ' n (p(t) -p{t j>n ))a(t) dW t + 

7=0 V-^j.n 
J r 2ir 



y+l,n 



a 2 (t)dt 



h N (t)a 2 (t)dt. 



Denote 



F N-=Y,J t {h N (t jtn )-h N (t))a 2 (t) 



dt. 



j=0 Jt i- n 



it holds 



< ^ 2sup|M*) - h{t)\+u h (p(n)) \E 



2 TV 



a 2 (t)dt 



2tt 



a 2 {t)dt 



Therefore by the condition p(n)N 2a — > oo as n — > oo we have 
(p(n))^/2 jB[ j j p 7V | M0 



as n — > oo. 



□ 
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Proof of Theorem 4.1. By the representation result in Theorem 4.2, 
we study the convergence in distribution of 



(p(n))-V 2 2 ^ h N (t j>n ) / ' / 



We write 



2 £ M*i,») / ' ( / dW s \a(t) dW t =: A 1 + A 2 



where 



-4, 



E (M*j>) - ^(<i,n)) [2 ( f o{s) dW s )a(t) dW t 

j—Q L Jtj,n ' 



and 



k " 1 r rtj+l.n / ft \ 

A 2 := £ 2 / ( / ff(s)dW s J 

Consider the term A 2 . Using Theorem 5.5 in [27] and generalization in 
Proposition 1 in [37], we have that 

(PW)" 1/2 E 2 / / <T(a)dWJ<7(t)dW t 

converges to 



2tt 



2H'(t)h(t)a 2 (t)dB t , 



where -B(t) is a Brownian motion independent from p(t) and the convergence 
is in law stably (for the definition of stable convergence see Section 2 in [27]). 
Therefore as p{n) — > the limit of 



k n 1 



(pW)" 1/2 E /l (^>) 2 / ' ( / a{s)dW s )o{t)dW t 

j—Q L ^ tj.n \Jtj,n / 

is mixing Gaussian with a random variance equal to 2 / 2?r H\t)h 2 {t)a A {t)dt 
and the Gaussian variable is independent from p{t). Finally consider the 
term A\. By (13) and Burkholder-Davis-Gundy inequality we have 



(P(n)r 1/2 E 



E ( h N(tj,n)-h(t hn ))( f tj+1 ' n ( f a(s)dW s )a(t)dWt 
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2 



< 



(pwrW sup \h N (t)-h(t)\ 

\0<t<2n 



X E 



E 

L \ j=0 - "i 



'j + l,n 



a(s)dW s )a(t)dWt 



<{ P {n)Y l (u h 



CE 



2- 



a 4 (t)dt 



where C is a suitable constant. Therefore (p(n)) 1 ' 2 A\ goes to in proba- 
bility as n — > oo. □ 

Proof of Theorem 4.5. The proof is similar to the univariate case, 
except it is necessary to consider nonsynchronous observations for different 
assets. With the notation in Theorem 3.3 we have to consider 



(40) 



M 1 

|fc|<iv v 



~N 



Ck(h)a k (N,pl,Pn) - E c k(h)c k (T, 

\k\<N 



12> 



We prove that 



(pwr 1/2 E i 



|fc|<7V 



N 



(41) 



(pH)~ 1/2 E (1 

|fc|<7V V 



c fc (/i)a fe (iV,p^,^) 
\k\ 



c k (h) 



^^D N {t)-t})f ti 



f 1 ,. At 2 , . 



t x Vi 2 , 



+ Op(l). 

In fact, we split 

1 1 



(42) 



2vr 2iV + 



T E E e^e-^-^Sji (p 1 )^ (p 2 ) J {/ i 



|s|<JV ij 

(43) + S E^^-^^^)^^^}- 

|s|<AT ij 
Consider (42). Observe that 

(44) J {/ i n</ 2 #0} = 1 {tlvt*<t} +1 /\t* +1 y 
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Therefore 



(45) 






(46) 




(47) 


+( 







1 At 2 



t 1 At 2 



1+1 dpHs)) (f % 3 dp\s)+ 3+1 dp 

J .7 1 + 1 .J + l 



t 1 , , Atf 



dp 2 & 



dp 1 (s) + 



■i+i 



t 1 At 2 
-' + 1 j + l 



We first consider (46). The terms (43) and (47) are estimated analogously. 
By splitting I| t i vt 2 <i i +iAt 2 j into the cases J/ + i W 2^ + 2 ^ + i 



{*i <* ? « ?+ 1 <* i+ 1 > ' {*i <*? <* J+i <* j+i > ' 



^{t 2 <t 1 <i 1 , ,<t 2 J and %?<*!<*=?. we obtain: 

e f 1 - 1) <*w e «^ Av(*? - 1, 1 ) ( d^ 

|k|^A NJ M U * 

dp\s)+ dp 2 {s)\l {t]yt]<t]+iMU} 

L j """i+l'^j+l 7 

2^ m<? w<? - 1\) dpHs)) (/; d^( 



8 + 



dp 2 (s) 



Define 



I7(0i,&) -g^Ar^Div^-^^i^j^i)^!^^ 
where Dpj(s) is the rescaled Dirichlet kernel. Then 

£M*?)Av(4-i?) f i+1 d P \s) f l dp 2 (s)I {t 2 <tl<tl t 

[/(01,0 2 )d P 1 (</» 1 )dp 2 (^ 2 ). 



2<<^1 

By ltd energy identity and Plancherel equality, we have 

2i 



2 1 
j+ll 



4>i<4>2 



U(<t> U <h)dp\fa)dp 2 {<l>2] 
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2- 



E 



<CiE 



D 2 N ((j)i - <fe)frjvW>i) £ %#2 



JO 



x ess sup ||£ 22 ||r/x> ess sup ||S 

2vi" 1 1 oo ii M 1 1 M 

- C2 2ivTT eSSSUP " S esssup ||E ||lo 
where C\ and C2 are suitable constants. It follows 

p(n)- l l 2 E 

< Cp{n)- X/2 N- X . 
Finally, we have to consider the term (45). Define 



Y,h N {t])D N {t}-t]) f ii+1 dp x {s) f l dp 2 (s)I {t2<t i <tl <t2 j 
• • Jt- Jt. 3 1 !+ J+ 



2t 



Vf" 



%i Ai +1 , 2 . . 

dp 2 (i). 



We have 



EMWffi)-!)/ 



Z 1 Af 
i+1 j+1 



t. 1 , 1 At j 



< p(n) AT esssup ||£ | |x, 00 esssup ||S ||l°c. 
It follows 

p(nr^E[\G N \ 2 ]<Cp(nf/ 2 N 2 , 
which goes to by the hypothesis. Therefore (40) is equal to 



EM*?)/. 



vt 2 

* 3 



dp\t) 



tWt 2 



dp\t) 



2tt 



h N (t)Y, l2 (t)dt + o p Q.). 



Then the proof proceeds as in Theorem 4.2. □ 
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